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Self-assembling novel ordered structures with nanoparticles has recently received much attention. 
Here we use computer simulations to study a two-dimensional model system characterized by a 
simple isotropic interaction that could be realized with building blocks on the nanoscale. We find 
that the particles arrange themselves into hexagonal superstructures of twin boundaries whose su- 
perlattice vector can be tuned reversibly by changing the temperature. Thermodynamic stability is 
confirmed by calculating the free energy with a combination of thermodynamic integration and the 
Frenkel-Ladd method. Different contributions to the free energy difference are discussed. 

PACS numbers: 61.50.Ah, 02.70.Ns, 64.70.Rh, 61.44.Br. 



The synthesis of nanoparticles and colloids with vari- 
ous types of interactions ilj and shapes [2] has advanced 
rapidly. This opens up the possibility to use them as 
building blocks for self-assembling novel materials 'from 
bottom up'. For future applications, ordered structures 
are especially interesting because they have unique pho- 
tonic and electronic properties. Furthermore, it is desir- 
able to have a large pool of different types of crystal struc- 
tures to choose from. Structure formation of nanoparti- 
cles is guided by the same principles as crystallization 
of atoms despite the different character of the interpar- 
ticle forces. Like for the majority of chemical elements, 
the most common ground states of nanoscopic building 
blocks are close-packed lattices which maximize density 
and the number of neighbor bonds. Different and more 
complex mesocrystals [3j have been observed experimen- 
tally for example in binary mixtures of colloids [H [5], 
cylindrical liquid crystalline phases [BJ, and in simula- 
tions of hard tetrahedra [7]. In general, however, the 
relationship between particle shape and interaction on 
the one hand, and the stabilized crystal structure on the 
other hand is not well understood. 

Among the most complex ordered phases are commen- 
surately and incommensurately modulated crystals [5], 
which have a basic structure and a superstructure in form 
of a spatial modulation. Although frequently found on 
the atomic scale in the bulk [9j [10] and as ordered struc- 
tures of noble gases in adsorbed layers [TTJ Q2], modu- 
lated phases have not been known to self-assemble with 
nanoparticles. So far only one-dimensional modulations 
have been forced upon a two-dimensional system of col- 
loids [13] and binary hard-disks [14] using an external 
periodic potential. The purpose of this work is two-fold: 
First, we show that planar superstructures with a tun- 
able modulation period form in a simple model system 
without external potential. Second, we confirm the ther- 
modynamic stability of the superstructures with free en- 
ergy calculations that are sensible enough to distinguish 
different modulations and analyze the origin of their en- 
tropic stabilization. 

Our model consists of identical particles interacting in 



two dimensions via a Lennard- Jones-Gauss potential |15] 
with parameters a 2 — 0.042, e = 1.8, r = 1.42. This po- 
tential is isotropic and consists of a single minimum with 
a shoulder (Fig.[]Ja)). Possible experimental realizations 
are spherical macromolecules and micelles. In the latter, 
a shoulder can be caused by a repulsive hard core of car- 
bon rings and an attractive soft alkyl corona [IB] E] ■ We 
integrate the particle motions with molecular dynamics 
(MD) computer simulations in the NPT ensemble at ex- 
ternal pressure P = 0. A potential cut-off at r = 2.5 
and a thermostat for temperature control are used. Typ- 
ical simulation times are 10 6 — 10 9 MD steps. Units are 
dimensionless and fixed by the choice of the potential. 
An exception is the temperature T, which is measured 
relative to a critical temperature Tq = 0.35. It has 
been shown in a previous study that the system trans- 
forms into a square crystal at Tq before it melts above 

T M = 0.40 ng. 

A system of 1024 particles is slowly cooled down using 
open boundary conditions to find the energetic ground 
state. Below 70% Tc, the system completely orders into a 
crystal with symmetry group p3ml (Fig.JTJc)). Since the 
particles are arranged in one triangle tile (Tr) and three 
slightly deformed pentagons tiles (Pe) per unit cell, we 
call this phase PeaTr. Additional squares (Sq) tiles (with 
two triangles attached) are present as defects. They ar- 
range along a line forming a twin boundary (TB) marked 
by an arrow in the figure. The effect of the TB on the 
crystal structure is a flip of the triangle tile orientations. 
We can measure the defect energy of a TB by comparing 
the potential energy of a perfect PeaTr crystal with one 
that has a single TB. The energy change per unit length 
caused by the presence of the TB is small compared to 
the average kinetic energy at the critical temperature, 
but positive: E/l = +0.13(1). This confirms that the 
creation of the TB requires energy. 

Next, we switch to a larger system with 10000 particles 
and study the temperature dependence of the equilibrium 
crystal structure in more detail. Fig. [2ja) shows the sys- 
tem held for a long time at 75% Tc- We used a rapid 
quench to T = to remove thermal noise. At the tern- 




FIG. 1: (color online) The interaction potential U(T) with 
a broad minimum (a) arranges the particles (disks) into 
(Pe)ntagons, (Sq)uares and (Tr)iangles (b). (c) The ener- 
getic ground state is PeaTr. The inset outlines a unit cell. A 
twin boundary is marked with an arrow. 

perature 75% Tc, TBs are found to appear in equilibrium 
and arrange into a highly ordered hexagonal superstruc- 
ture (symmetry group p6m) : (i) The orientation of TBs 
is restricted to three directions, (ii) TBs are in average 
equally spaced, (iii) TBs form triple intersections only, 
double intersections never occur. These observations are 
sufficient to identify the superstructure as an interface 
modulated |19j or, alternatively, domain wall incommen- 
surate phase [i"2] . 

Given the simple pair interaction, the superstructure 
in Fig. [5Ja) is surprisingly complex. The average period- 
icity equals 33 nearest neighbor distances and the super- 
structure unit cell contains approximately 700 particles. 
Furthermore, as will be shown below, this superlattice 
constant can be tuned continuously by varying the tem- 
perature. For this to happen there needs to be a way for 
the TBs to propagate efficiently through the structure in 
form of a phason mode [5] . The basic mechanism for such 
a motion is depicted in Fig. [2jb) . Squares can change 
position in three steps: (i) The combination of a square 
and a neighboring triangle transforms into a pentagon by 
moving the two particles of the common edge outwards. 
As a result, an isolated triangle remains, which shares 
edges with three pentagons, (ii) This transformation is 
energetically disfavorable and generates internal stress, 
which can be relaxed by transforming another pentagon 
nearby into a combination of a square and a triangle, 
(iii) Subsequent propagation of all squares in a TB in- 
duces a motion of this TB. 

The hexagonal superstructure is best analyzed in 
Fourier space. Diffraction images (DIs) at six differ- 
ent temperatures are shown in Fig. [3] Up to 70% Tc 
the PesTr crystal is stable. The DI at this tempera- 
ture consists of equally spaced, sharp peaks. Only weak 
diffuse scattering caused by thermal motion is present 
in the background. Above 70% Tc, the peaks split up 
into a strong central peak and weak satellite peaks as 
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FIG. 2: (color online) (a) Ordering of twin boundaries into a 
hexagonal superstructure. Only the tiling, not the particles 
are shown, (b) A structural rearrangement can transform a 
pentagon into a square plus triangle and induce a motion of 
twin boundaries. This collective motion is a phason mode. 

expected for a modulated phase. The wave vector sep- 
arating the satellite peaks continuously increases with 
increasing temperature as a result of a decreasing super- 
lattice constant. A change in the rotational symmetry 
of the DI is visible in Fig. [3ja) : The six- fold symmetry 
present at low and intermediate temperatures gradually 
transforms into a twelve-fold symmetry (generalized sym- 
metry group p\2m). Note that at 97% Tc, twelve peaks 
form two dodecagonal rings in the DI as outlined by a 
dotted circle. This suggests that a dodecagonal qua- 
sicrystal is stabilized close to the critical temperature. 
The structure and dynamics of this quasicrystal has been 
studied recently |18j . The sequence of symmetry changes 
is summarized as: 

p3ml ' 7Tc > p6m °' 97Tci > pYlm p4m — ^> E(2). 

(1) 

While the first two phase transitions are second order, 
the transition to the square crystal at Tc and the melting 
transition at Tm = 1.15Tc are first order [18] , Snapshots 
of the particle configurations at various temperatures and 
full diffraction patterns can be found in the supplemen- 
tary information. 

We determine the superlattice wave vector by measur- 
ing the distance q of two neighboring superlattice peaks. 
Simulations are started from PesTr (stable at low T) and 
the fluid (stable at high T) to test for hysteresis in the 
temperature behavior. The dependence of the superlat- 
tice wave vector on temperature is plotted in Fig. [4] No 
hysteresis is present and the final configuration is found 
to be independent of the starting condition, which is 
an indication that thermodynamic equilibrium has been 
reached. The wave vector of the superstructure increases 
roughly linearly between T = 0.24 « 70% T c and T c . 
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FIG. 3: (color online) (a) Each sector shows a diffraction 
image obtained by fast Fourier transform of the time aver- 
aged particle density. The temperature increases counter- 
clockwise, (b) Magnified views of slices of the diffraction im- 
ages (indicated by a dashed rectangle in (a)) show the increase 
of the superlattice wave vector. 

Such a continuous increase of the wave vector is typical 
for two-dimensional systems and is characteristic for a 
'floating' modulation [S]. 

The stabilization of the superstructures is analyzed in 
more detail by calculating free energies. Two comple- 
mentary techniques are used: the Frenkel-Ladd method 
(FL) and thermodynamic integration (TI). The absolute 
free energy can be obtained at all temperatures combin- 
ing the methods into a new method (FLTI): (i) At low 
temperature (< 10% Tc), the absolute free energy is de- 
termined with FL by interpolating between the target 
structure and an Einstein crystal with known free en- 
ergy • (ii) The free energy is then extended to higher 
temperatures with TI by integrating the potential energy 
E(T) along a reversible path: 

F(T) = c(T )T-T£dT>^p_, (2) 

where To is a reference temperature and c(Tq) an inte- 
gration constant that is determined from a comparison 
with FL. Using both methods is necessary because FL 
cannot be applied at high temperatures. The reason is 
the motion of the TBs, which makes equilibration during 
the interpolation to the Einstein crystal very slow due to 
the presence of particle rearrangements. 

Boundary conditions are important for equilibrating 
the superstructures. Open boundaries ensure that TBs 
can be created and annihilated. In contrast, by switch- 
ing to periodic boundary conditions, the superlattice con- 
stant q remains fixed because the total number of TBs in 
the system cannot change. Using periodic boundaries al- 
lows to simulate superstructures at temperatures where 
they would otherwise be metastable with regard to a 
change in q. 




Temperature T 

FIG. 4: (color online) The wave vector of the superstructure 
grows linearly with temperature. No hysteresis is observed in 
molecular dynamics simulations. The growth of the superlat- 
tice wave vector is confirmed by free energy calculations. 

We analytically construct superstructures with various 
values of q and calculate the free energy F(T,q) using 
FLTI and periodic boundaries. Stable superstructures 
are characterized by dF/dq\x = 0. Results from free 
energy calculation are included in Fig. [4j They agree 
well with the MD simulations and confirm that the tem- 
perature dependence of the superstructures wave vector 
observed in MD simulations. 

The stabilization of the superstructures at elevated 
temperatures can only be explained by a higher entropy 
in comparison to the PesTr crystal. We now analyze the 
contributions to the free energy. In general, the free en- 
ergy of a classical <i-dimensional system can be written 
as a series: 

F(T) = E(0) - §Nk B T\n(T) + c'T + aT 2 + (3T 3 + . . . 
. ' . ' 

F H (T) F A (T) 

(3) 

with new constants c', a, and j3. This series splits up 
into a 'harmonic' part Fn(T) and an 'anharmonic' part 
Fa(T). The harmonic part corresponds to the free energy 
of the system in a linear approximation, where the dy- 
namics is completely specified by the dynamical matrix, 
or alternatively, by the zero temperature phonon density 
of state. The constant c', and therefore Fh, can be de- 
termined from the FL results at low temperature. The 
anharmonic part contains higher-order effects and has a 
more complicated temperature dependence. We write Fa 
with the help of Eq. ([2| in a form that can directly be 
used in simulation: 

Fa(T) = -T J dT'^^ \^ (4) 

Note that since the integral does not diverge for small 
temperatures, we can now take T = as the lower limit. 

In order to study the temperature dependence of Fa., 
we plot the integrand of this equation for the low temper- 
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simulation techniques and theoretical ideas outlined in 
this work are of interest. The combination of the Frenkel- 
Ladd method and thermodynamic integration is a pow- 
erful tool that allows to calculate free energies with very 
high precision (< 10~ 4 ). In our case it was possible to 
compare superstructures with different wave vectors and 
distinguish the free energy contributions from oscillatory 
(phonon) particle motion and phason modes. 

Support from the Deutsche Forschungsgemeinschaft 
(EN 905/1-1) is gratefully acknowledged. 



FIG. 5: (color online) Analysis of the entropic stabilization 
of the quasicrystal and the superstructures. The anharmonic 
part of the free energy per particle, Fa(T)/N, can be obtained 
by integrating [E(T) - E(0) - Nk B T]/NT 2 , plotted in the 
figure. The dashed line splits contributions originating from 
oscillatory (phonon) motion and the phason modes. 

ature phase PesTr (CR) and the quasicrystal (QC) stable 
at high temperature in Fig. [5j The shaded area between 
the curves corresponds to a free energy difference of the 
two phases. 

The free energy difference at the critical temperature, 

AF = [F QC (T C ) - F QC (0)] - [F CR (T C ) - F CR (0)], (5) 

has three contributions: (1) Harmonic particle motion 
in the local potential energy minima (HM). (2) Anhar- 
monic particle motion in the local potential energy min- 
ima (AM). These two parts vary smoothly with temper- 
ature. (3) Particle rearrangements as part of the pha- 
son modes (PH), which correspond to the TB motion, 
but can also be pentagon rotations [T5]. While the first 
contribution is equal to TpAc', the last two are calcu- 
lated from the shaded area in Fig. [5| Notice that pha- 
son modes only contribute after they are activated at 
around 70% T M . We estimate their contribution by ex- 
trapolating the low-temperature data for the quasicrys- 
tal as indicated by the dashed curve. The results are: 
AF RM /N = 0.042, AF AM /iV = 0.017, AF PH /N = 0.010. 
The majority of the free energy difference comes from os- 
cillatory (phonon) motion. Interestingly, phason modes 
do not contribute significantly. 

In conclusion, we have introduced a model system 
that stabilizes a continuum of modulated superstruc- 
tures. The interaction potential is simple enough that 
it could in principle be realized for example with attrac- 
tive nanoparticles. Complex two-dimensional structures 
are now feasible experimentally: Only recently complex 
tunable two-dimensional binary molecular networks have 
been fabricated [2Tj . In another work, it has been demon- 
strated that binary nanoparticles can be assembled at the 
liquid-air interface into highly ordered crystals [22 . How- 
ever, even without direct experimental realization, the 
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Entropic Stabilization of 
Tunable Planar Modulated Superstructures 
- Supplementary Figures - 



Figure SI. Modulated superstructure at temperature T= 0.25, close to the temperature 
where the superlattice constant diverges. Despite the large distance of the twin 
boundaries formed by squares (magenta), they form a well-defined hexagonal lattice 
with a huge unit cell (several thousand atoms). Note that the twin boundaries cross in 
triple intersections only. 
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Figure S2. At temperature T = 0.27, the superlattice constant has decreased compared 
to T = 0.25. This is possible by the formation of new twin boundaries at the surface 
(open boundaries) and their motion into the system via a flip process. The 
crystallographic orientation of neighbouring single crystalline domains separated by a 
twin boundary is related by a mirror operation along the twin boundaries. 
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Figure S3. The dodecagonal quasicrystal grown at T= 0.30 and then quenched to T— 0. 
Due to the periodic boundary conditions used in the simulation, new twin boundaries 
cannot form and the system cannot transform into a superstructure. It is observed that 
squares are always accompanied by two triangles. 
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Figure S4. Diffraction patterns at various temperatures (part 1). From top left to bottom 
right: T= 0.20, 0.25, 0.27, 0.29. 
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(c) <K> ~ <X> ~ O 

Figure S6. The tiling of the superstructures and the quasicrystal can be mapped on 
another tiling by connecting the centres of the pentagons (cyan). The resulting new tiles 
are triangles, squares, shields, and an acute ellipsoidal tile obtained by the intersection 
of two dodecagons. The shield is centred by a single triangle (yellow) of the original 
tiling; the ellipsoidal tile by two triangles and a square (magenta), (a) The alternative 
tiling around the intersection of three twin boundaries, (b) The alternative tiling for the 
section of the dodecagonal quasicrystal. (c) A local rearrangement in the alternative 
tiling. 
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Figure S7. Free energies of superstructures with various superlattice wave vectors (Aq) 
as a function of temperature. The intersection of two curves for different superlattice 
wave vectors is used to estimate the equilibrium transition temperature between 
superstructures with these two wave vectors. The two lines correspond to calculations of 
the free energies in a heating run (starting from the Pe3Tr crystal) and a cooling run 
(starting from the liquid). 



7/7 



